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ABSTRACT 

The two main characteristics of the holmium ions in H0F3 are that their local 
electronic properties are dominated by two singlet states lying well below the 
remaining 4f-levels, and that the classical dipole-coupling is an order of magni- 
tude larger than any other two-ion interactions between the Ho-moments. This 
combination makes the system particularly suitable for testing refinements of 
the mean-field theory. There are four Ho-ions per unit cell and the hyperfine 
coupled electronic and nuclear moments on the Ho-ions order in a ferrimag- 
netic structure at Tq = 0.53 K. The corrections to the mean-field behavior of 
holmium triflouride, both in the paramagnetic and ferrimagnetic phase, have 
been calculated to first order in the high-density l/z-expansion. The effective 
medium theory, which includes the effects of the single-site fluctuations, leads 
to a substantially improved description of the magnetic properties of H0F3, in 
comparison with that based on the mean-field approximation. 



PACS-numbers: 75.10.-b ; 75.30.-m ; 75.50. Gg 



I. Introduction 

The magnetic properties of H0F3 have been established experimentally 1-3 in 
considerable detail, and these experiments were recently interpreted in terms of 
a mean- field (MF) model by Leask et al. 1 - H0F3 is orthorhombic, and the low 
symmetry at the Ho sites implies that the J = 8 ground-state multiplet splits 
into singlets. In the paramagnetic phase, the energy difference between the two 
lowest singlets is about 0.7meV and the dipole matrix-element between the two 
states is large. The next singlet lies at about 5 meV whereas the remaining ones 
lie between 10 and 50 meV above the ground state, 4 and this arrangement of the 
crystal-field levels leads to a very anisotropic susceptibility at low temperatures. 
The system is close to an Ising system, with the only modification that there 
are two easy directions, one for each magnetically-equivalent pair of Ho-ions 
in a unit cell. The classical dipole-coupling between the angular moments is 
weak, but is nevertheless found to be one order of magnitude stronger than 
any other interaction between the dipoles. It gives rise to strong correlation 
effects below 5K and is responsible for the induced magnetic ordering of the 
singlet ground-state system at Tq = 0.53 K. The dipole coupling is however 
not sufficiently strong to produce an ordering of the isolated electronic system. 
The ratio between the two-ion coupling and the threshold value required for 
inducing magnetic ordering of the electronic moments is found to be R = 0.86. 
The phase transition occurs only because the hyperfine interaction between 
the electronic and nuclear moments enhances the effective susceptibility, thus 
leading to a cooperative ordering of the two systems. Below Tq the ordered 
moments are along the two easy directions. At low temperatures, H0F3 may be 
considered to be a simple singlet-singlet system in which the moments interact 
like classical dipoles. The dipole-dipole interaction can be calculated accurately 
from first principles and the magnetic system is fully characterized by only a few 
parameters. The only complications are that there are four magnetic ions per 
unit cell, that the hyperfine interaction plays an active role, and that additional 
two-ion couplings are of some importance. 

Although the mean-field theory of Leask et al} was able to reproduce many 
of the observations made in H0F3, such as the excitation spectrum determined 
by neutron scattering at 1.6 K, some discrepancies remained. The calculated 
moment in the zero temperature limit was 16% smaller than that observed, 
and the comparison between the theoretical and experimental heat capacity 
was not entirely acceptable. Quadrupole-quadrupole couplings may be of some 
importance in H0F3 for explaining the discrepancies, but in this paper we shall 
concentrate on the corrections to the mean-field theory due to the single-site 
fluctuations. In section II we recapitulate the results of the high-density 1/z- 
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expansion for the singlet-singlet model utilizing the effective medium approach. 5 
The theory is then extended to cover both the paramagnetic and the ordered 
case. In section III the theory is applied to H0F3 and the conclusions are 
presented in the final section IV. 

II. The 1/z-expansion 

We consider a Bravais lattice of N identical singlet-singlet atoms, characterized 
by the energy separation A between the two singlets, and by the dipole matrix- 
elements. In the case of a singlet-singlet system, the x-axis may be defined so 
that only the matrix-elements of J x are non-zero, and M denotes the (numerical) 
value of this J x matrix-element between the two states, whereas +m and — m 
are the matrix-elements of J x within respectively the lower and upper singlet. 
Mq = M 2 + m 2 is a constant and m = in the paramagnetic case at zero field. 
At the temperature T = l/k B (3, the MF population-factors of the lower and 
the upper state are n and n 1 , where n = 1/[1 + exp(— j3A)}, n + n x = 1 and 
we define n 01 = n — n l . 

The two-site Green function is defined as the r-ordered ensemble average 

G(u,ri - r 2 ) = -(^1^)1^)) (2.1) 

\ti — r 2 | < (3, and we use the short-hand notation 

Jix = Jix ~ ^x) (2-2) 

The Fourier transform of the Green function is defined in terms of (h times) the 
Matsubara frequencies, u n = 2nn/j3, where n is an integer: 

G(q, iw n ) = E e- tq ' R - f G(ij, r) dr (2.3) 
. JO 

The Hamiltonian is divided into two parts, H = Hq + Hi, where TCq is the 
mean-field part and H 1 is the perturbation 

«i = -5E^')4^ (2-4) 

y 

in which case the Green function is determined by the linked-cluster expansion 



with 



U(P,0) = 1- /^i(r 1 )rfr 1 + --- (2.6) 

J 



The index on the thermal averages indicates that they are the mean-field 
values. The non-interacting Green function, obtained in the zeroth order of Ti-^, 
vanishes if i ^ j, whereas if i — j it is 



o(^n) = -M 2 g(iu n ) ~ m 2 h{iuj n ) 



(2.7) 



where the two response functions are 




h{iu n ) = (3(1 -ng^^o 



(2.8) 



A 2 - (iu n ) 
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The perturbation Tii cannot, in general, be considered as being small com- 
pared to 7i , but each time a term involving the two-ion coupling is summed 
over q, we effectively gain a factor 1/z, where z is the co-ordination number. 
A systematic expansion of the Green function in powers of 1/z, in the case of 
the singlet-singlet system, was developed by Stinchcombe. 6 Here, we shall use a 
slightly different approach and utilize the concept of an effective medium, which 
is the basis of the coherent-potential approximation. 

The ensemble averages in (2.5) are calculated by expressing the angular mo- 
mentum operators in terms of the 'standard-basis' operators, a vp . = \v){fJ>\ where 
\u) are the MF-eigenstates of an atom, and by utilizing the invariance of the 
trace to a cyclic permutation of the operators. 7,8 These are not Bose-operators, 
so the 'contractions' determined by the commutators of the different operators 
are not c-numbers, but operators which give rise to new contractions. However, 
the next generation of contractions adds terms to C(q, iu n ) which always involve 
additional q-summations. Hence in the order (l/z)°, these contractions are ne- 
glected. This corresponds to a decoupling of the higher-order cumulants in (2.5) 
into products of the second-order terms, (T T J \ X (T]) J j x ( T 2)) o = ~~ ^j^o( r i ~~ T 2) 
(if the possible difference between (J x ) and (J x ) is neglected, see below). The 
infinite series of 'chain diagrams', generated by (2.5) using this decoupling, is 
easily summed, and the result is: 



showing that this approximation is equivalent to the random-phase approxima- 
tion. The difference between the fourth-order cumulant and the corresponding 
decoupled value appears in the next order of 1/z. In the usual 'unconditional' 
cumulant expansion this difference, the fourth-order semi-invariant, is intro- 
duced as an additional vertex. The vertex in the RPA chain-diagrams is re- 
placed by the sum of this and the fourth-order semi-invariant, and neglecting 
any particular effects of whether some vertices in a chain belong to the same 
site or not, we may straightforwardly sum the series. 6 





(2.9) 



RPA l + J(q)G (iu n ) 
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In the order 1/z only the single-site Green function is directly modified, as 
the fourth-order cumulant only differs from the decoupled one if all the four 
momentum operators belong to the same site. The single-site Green function is 

G(iu> n ) = G(ii,iu> n ) = ^£G(q,^J (2.10) 

and, to the order 1/z, the two-site Green function may be expressed in terms 
of G(iu> n ) by introducing the effective-medium coupling 

= ^E^(q)G(q,^ n )/G(^J (2.11) 
JV q 

in which case 

G{ito ) 

*^ = i tW -™ (2 ' 12) 

K(icu n ) is the sum of all chain diagrams which start and end at the same site 
without crossing this site in between. 5 We may therefore consider an effective 
cumulant expansion of G(iu> n ), equivalent to (2.5) where J {if) in Ti 1 is replaced 
by the time-dependent coupling K{ji — r 2 ), and the term to leading order in 
this coupling is 



-iff I 

A Jo Jo Jo 



P „i 



X 



( T TJix( T 2)Jix(T 1 )J lx ( r )Jix(Q))o - ( T r^( r 2)^(^l))o( T r^( r )^(0))o 

(2.13) 

The decoupling of the higher-order cumulants in the single-site series into prod- 
ucts of the second-order terms leads to the result 

°^U = rnSfcj (2 ' 14) 

and when this result is introduced in (2.12), K{iu n ) cancels out and we get the 
RPA-result, (2.9). To the next order in 1/z we have to include the fourth-order 
cumulant term, (2.13). We shall first consider the paramagnetic phase, m — 0, 
in which case 

G(iu n ) = G Q {iu n ) - G (iu n ){K(iu n )G (iu n ) + S(i^ n )} + • • • (2.15) 
where the renormalization factor is 

= a + ~f{iu n )g(iu n ) (2.16) 

a is a constant 



a 



n oi 



*2-\{9{G) + P{l-r$i)}\i\ (2.17) 
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and the frequency- dependent term is 



M 2 r i 

7 (^J = -5- Ai - (1 -n 2 01 )K(iu; n )] (2.18) 



n, 



01 

and the parameters \ p are defined as 

\ = -oY. K ^n')W^)] P (2-19) 
" n' 

The unconditional cumulant expansion accounts correctly for the fourth-order 
cumulant term in (2.15), but an analysis of the sixth and higher-order terms 
shows that this procedure does not lead to a good estimate of the higher-order 
contributions in the single-site series. Instead, it is found that the series gen- 
erated by replacing G (iu n ) in front of the second term of the single-site series 
(2.15) by the interacting Green function G(icu n ), much more effectively accounts 
for the terms deriving from, specifically, the sixth-order cumulant. Introducing 
this Dyson-like result for the single-site Green function in (2.12), we get: 

a ^ = TTW^m^3 (2 ' 20) 

valid to first order in 1/z. This result is nearly the same as the one derived 
by Galili and Zevin 9 using an elaborate renormalization of the unconditional 
expansion. In addition to the simplifications attained by utilizing the effective- 
medium approach, the present procedure is fully self-consistent. A more detailed 
discussion of the effective medium theory and its comparison with the uncon- 
ditional cumulant expansion may be found in Refs. 5 and 10. In a systematic 
expansion in 1/z, the effective-medium approximation ceases to be valid in sec- 
ond order. However, an improvement of the theory is obtainable by including 
the new diagrams due to the sixth-order cumulants into the single-site series, 
neglecting the corresponding {1/z) 2 two-site effects. The second-order contri- 
butions to the effective medium have been considered in an analysis 11 of Pr and 
are of some importance if the energy gap in the excitation spectrum is small. In 
the present analysis we shall consider only the leading-order 1/z- modifications 
due to the single-site fluctuations. 

In the paramagnetic phase, the contribution of the fluctuations to the in- 
ternal energy is 

W = ^ E {K(iu n ) - ^p±}G{iu n ) - Nn,A (2.21) 

and introducing the l/z-result for G(iu> n ) in this expression, we may write the 
energy change: 



N 



a7l 01 A a ,2- 1 



A - M% + - £ K(iu n ){G(iu n ) - G (iu; n )} (2.22) 



1 + a 
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where the last sum is a small second-order term. The derivative of 5U with 
respect to T determines the change of the heat capacity from its mean-field 
value. In the calculation of the energy it is of importance that the Green 
function satisfies the sum rule: 

WG^n) = -M 2 (2.23) 

P n 

(m — 0). This is strictly the case if the single-site series is terminated after the 
first or the second term in (2.15). It is also found numerically to be valid with 
a rather high accuracy for the self-consistent 1 / z-Green function. 

Next we wish to consider the situation when m^fl. In the low-temperature 
regime, with which we shall mostly be concerned, the elastic pole in (2.7) is 
weak compared to the inelastic one, unless the system is close to the saturation 
limit m 3> M. However, in this limit, i.e. at high fields, all the fluctuations 
are quenched and the mean- field approximation is valid. This means that, 
generally, at low temperature the elastic fluctuations are of less importance 
than the inelastic ones, and it is therefore acceptable to include the elastic-pole 
contributions less rigorously. 

The first problem we meet in the order 1/z is that, if we define the mean- 
field Hamiltonian in the usual way, (J x )q will in general be non-zero, giving 
rise to additional complications in the linked-cluster expansion of the Green 
function. In order to avoid these complications we introduce a modified mean 
field, H MF , defined by the requirement that 

(J x ) = (J x )o = m 2 n 01 (2.24) 

which differs from the usual mean field H : 

H MF = H -5H; gv B H = gn B H + J(0)(J x ) (2.25) 

where H is the applied field. The determination of Hq in terms of -f^F instead 
of H Q introduces an extra perturbation in the Hamiltonian, Ti = Hq + H 1 +H 2 , 
with H 2 = ~9l x B&H*52iJ xi- The r& tio 5H/H MF is a small quantity, at the 
most about 0.1 in the numerical calculations, and 7Y 2 is truly a weak perturba- 
tion. Furthermore, the leading order modification of the cumulant expansion is 
proportional to (5H/H MF ) 2 , which small correction is neglected. 

The elastic pole disturbs the inelastic fluctuations, and H(iu> n ) in (2.20), 
when n^O, is changed into 

r 2m 2 i 

E(iu n ) =a-a m + [y(iu n ) - j^li )] 9 fan) ( 2 - 26 ) 

with 

m 



,2_ 4 



10, 



01 



A 2 - 9(0)\i + ^ A 3 - (1 - r&Xl + ±PAn 01 )K(0)g(0)\ (2.27) 
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In this result we have for simplicity neglected some frequency-dependent terms 
which vanish in the zero-frequency limit (g(iou , ± iuj n ) appearing as a factor 
in the n'-sums has been replaced by its zero-frequency value g(iu n ,)). The 
inelastic broadening of the elastic peak discussed in Ref.ll is neglected in this 
approximation. In addition to the inelastic modifications we get elastic terms 
appearing only at zero frequency, and at this frequency the single-site Green 
function is found to be 

G < )-- l + S (0) M M°i-(0) g (0) -'^'< » < 2 ' 28 > 

where E(0) is the zero-frequency value of Yi(iu n ) given by (2.26) and 

l+tanMm^^ffl-M 2 ^) 
5 l + {4nl 1 K(0)g(0) + 2\ 2 + g(0)\ 1 }Mynt 1 { ' ] 

The cumulant expansion of the elastic contribution effectively becomes an ex- 
pansion in (3K{Qi) which converges slowly at low temperatures, where however 
the elastic term as a whole is frozen out. The actual result for the 1/z-term 
in the numerator of £, deriving from the fourth-order cumulant, is the lead- 
ing order term in the Taylor expansion of the tanh-term, which diverges in the 
zero-temperature limit. This divergence was not found to influence the nu- 
merical results significantly, as the term is multiplied by h(0) which vanishes 
exponentially in this limit. However, in order to account in a simple way for 
the renormalization of this term expected due to higher-order contributions, we 
have replaced the 1 / z-term by its tanh- value. 

The final result for G(q, iuj n ) to first order in 1/z when m is non-zero is 
determined in terms of G(iu n ). At non-zero frequencies the result may be 
written as (2.20), with Y,(iui n ) now given by (2.26). At zero frequency we have 
to introduce G(0) given by (2.28) in the original expression (2.12) for G(q, iuJ n ). 
The result (2.12) becomes equivalent to the RPA-result if G (iu n ) in (2.9) is 
replaced by the effective non-interacting Green function 

G (iu n ) = ^%rr. r (2.30) 

which is equal to G (icc; n ) / '{1 + ^(^ u; n)}> except in the zero-frequency case when 
m ^ 0. 

The mean field H MF is determined by (2.24) which may be written 

(4) = Wo = -9VB / G (0; H')dH' (2.31) 

where G (0;H') is the non-interacting Green function, (2.7), as a function of 
the Zeeman-field H' . Differentiating this equation with respect to the applied 
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field g/2 B H, at a constant temperature, we get on the left-hand side the static 
bulk-susceptibility —(2(0,0), or 

G(0, 0; H) = G (0; H MF ){1 - J(0)G(0, 0; H)} (2.32) 

dH 

which is directly integrable, as H MF is the mean field determining G(0,0;if), 
and we get the following relation between the two field quantities 

f H MF G n (0;H') 

H Q = / ! -dH' 2.33 

Jo G (0;H>) 

making use of the effective non-interacting Green function defined by (2.30). 
The combination of (2.33) and (2.25) determines the relation between the ap- 
plied field and the mean field, which allows a fully self-consistent calculation of 
(J x ) as a function of field, in much the same way as in the MF- approximation. 
The adjustment of the mean field, by replacing H by the effective value H^p, 
implies that the change in the free energy including the l/z-contributions, rel- 
atively to the non-magnetic state, is determined by the mean-field part of the 
Hamiltonian, 7i , because dF/N = —gfi B (J x )dH is equal to —gfj JB {J x ) dH (at 
constant temperature). This has the consequence that 

r M o 

5F(m = 0) = -N / gfi B 5Hd(J x ) (2.34) 

J 

Mq is the saturation value of (J x ) in the limit of an infinite field, in which 
limit any correction to the MF- approximation vanishes. The difference 5H = 
H — H MF is considered to be a function of (J x ), and introducing 5F(m = 0), 
as determined by 5U in (2.22) and the corresponding change in heat capacity, 
the fulfillment of (2.34) provides an independent test of the theory. 

III. The magnetic properties of HoF 3 

H0F3 is orthorhombic with the lattice parameters a = 6.404 A" 1 , b = 6.875 A -1 
and c = 4.379 A -1 , and there are 4 Ho 3+ and 12 F - ions per unit cell. The 
positions of the ions within the unit cell are specified in the previous papers 2 ' 3 
on H0F3, and the projections of the ions on the a-c and a-b planes are shown in 
Fig.l. The crystal-field Hamiltonian for the Ho-ions at the four different sites 
in the unit cell is the same when referring to one local coordinate system for 
the sublattices labeled 1 and 2 in Fig.l and another, in which the b- and c-axes 
are reversed, for the sublattices 3 and 4. The local easy axes, or rr-axes, are 
indicated in the figure. They lie in the a-c plane and make the angles 6 and 
—9 with the a-axis. The magnitude of 9, but not its sign, is known from the 
experiments. 1 ' 2 The sign of 9 has not much influence on the analysis, but we 
have here chosen 9 = +25° as also used in Fig.l. This choice of sign leads to a 
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slightly better fit of the magnetization curves than the other, but the difference 
is not so significant that it rules out the alternative choice. 

The equivalence of the ions, when considered in their local coordinate sys- 
tems, makes it straightforward to generalize the theory in the preceding section 
to this case with four magnetic ions per unit cell. In the paramagnetic phase 
the only modification of the MF/RPA-calculations presented in Ref.l are that 
the components in the single-ion susceptibility tensors x°(u;), deriving from the 
singlet-singlet transition on the rth ion, should be divided by the common renor- 
malization factor 1 + E(u;). This follows from the fact that the RPA-expression 
for the final Green function is unchanged if the non-interacting Green function 
G (iu n ) in (2.9) is replaced by the effective G (iu n ) given by (2.30). Neglecting 
small 1/z-corrections of the RPA-contributions from other single- ion transitions, 
this procedure allows the full level scheme of the ions to be included in the cal- 
culations. The imaginary part of £(u;) is proportional to the imaginary part of 
K(uS) times (1 — n^), and the latter factor is vanishing small at low tempera- 
tures, thus E(u;) is real just above Tq. In this limit the paramagnetic density 
of states, for each of the 4iV magnetic ions, is 

^) = T^E^q,,-^) (3-1) 

where e q v are the four excitation energies at a certain q. The q-summation in 
the expression (2.11) for K(iu n ) may be replaced by an energy integration by 
the introduction of N{e). The resulting integral is the same as the one derived 
for a Bravais-lattice. Calculating K{iu) n ) in this way, the final density of states 
has been determined self-consistently by an iterative procedure, using the RPA 
as a starting point. The result, nominally at 0.55 K, is shown in Fig.2. The 
corresponding excitation energies at 1.6 K along the high-symmetry directions 
(hOO) and (0 0/) are shown in Fig. 3. As discussed in Ref.l an effective double- 
zone representation is applicable along these directions, and only two modes 
scatter the neutrons in a constant-q scan. The calculated dispersion relations, 
the solid lines in Fig.3, are compared with the experimental neutron-scattering 
results and the RPA-predictions presented in Ref.l. Although the excitation 
energies in the order \jz are not changed much compared to the RPA-theory, 
the tendency of the excitations to split into two separated bands is much more 
pronounced in the final density of states in Fig.2 than in the RPA. 

Some of the parameters in the present calculations differ from those used in 
the previous MF/RPA-model. Including 8U given by (2.21) in the calculation 
of the paramagnetic heat capacity, the fit to the experimental position of the 
maximum in the Schottky anomaly is improved by a slight reduction of the 
zero-field value of A from 0.71 meV to 0.69 meV. This change is accomplish by 
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using 



V m = -3.42 10~ 5 meV 

V 62 = (6.41 -i 4.59) l(T 5 meV 



(3.2) 



whereas the remaining crystal-field parameters are the same as in Table I in 
Ref.l. The change also influences weakly the matrix element (M = 6.551 instead 
of 6.541). The quality of the fit to the experimental heat capacity is much 
improved in comparison with the one derived in the MF-approximation, as 
shown in Fig. 4. The fit is now close to perfect below 7K, where the phonons 
are unimportant, and down to about IK, where critical fluctuations may start 
to be of importance. 

The dominant part of the two-ion coupling, due to the classical dipole- 
dipole interaction, is unchanged. Given the positions of the magnetic Ho-ions, 
the dipole-dipole interaction is calculated by the Ewald method utilizing the 
refinements developed by Bowden and Clark. 12 In order to get agreement with 
the observed transition temperature, and to explain the larger splitting between 
the A\ and A 2 excitation modes at (10 0) than that predicted by the classical 
coupling, two exchange-coupling constants, and J^ii were used as fitting 
parameters. They are the isotropic couplings between respectively nearest and 
next-nearest neighbors. To these we here add one more coupling, J~ n , which is 
the nearest-neighbor coupling between ions on the same sublattice (the distance 
between these ions is not much greater than the smallest one). The final values 
of the three parameters are: 



where the numbers in the brackets are the previous values used in the MF/RPA- 
model (when 9 = +25°). The introduction of the coupling parameter J7q leads 
here to a slight improvement of the fit to the excitation energies, and thus 
probably also to a better estimate of the density of states, but it is in no way 
essential for the analysis. For instance, the reduction of the calculated energy 
of the upper A r mode close to (1 0) is a pure renormalization effect. The new 
values of the exchange parameters mean that the effective coupling parameter 
determining Tc (see Ref.l) 



is now 7.92 /xeV, which is about 11% larger than in the MF-model. This increase 
is required in order to compensate for the reduction, by the factor 1 + S(0), of 
the effective single-site susceptibility deriving from the fluctuations. This factor 
is calculated to be 1.128 just above Tc- 



3\z 
Ju 



0.30 /xeV (-0.64 //eV) 
0.37 /ieV (0.30 //eV) 
0.12/ieV (0) 



(3.3) 



J X {Q) = 2J 13 cos 29 + 2J 12 + 2J n + 16.848 J, 



(3.4) 
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In the calculation of the zero-frequency susceptibility, it is important to 
include the nuclear contribution due to the hyperfine interaction. This cou- 
pling enhances the susceptibility at Tq by about 16%. To a first approximation 
the influence of the fluctuations on this coupling may be neglected. However, 
it is practicable to include the effects of the singlet-singlet fluctuations on the 
hyperfine coupling, leaving out the intrinsic 1 / ^-modification of the nuclear sus- 
ceptibility (which is of the order A 4 ). In this approximation the rcrc-component 
of the effective no n- interacting susceptibility, x° — X°/{1 + ^(0)}> 1S replaced 
by 

XeS = X°(l + A 2 x°Xi) (3.5) 
and K(0) by approximately K(0)(1 + A 2 x°Xi)- I — 7/2 is the nuclear spin, 
A = 3.36 yueV is the hyperfine coupling and Xi is the nuclear xx-susceptibility 
component as determined from 7i hf (eff) given by (3) in Ref.l. x i n this 
Hamiltonian is replaced by j^SqXlq)) but this change can be neglected in 
the quadrupole term. 

The maximum value of the renormalization factor, 1 + E(0), is 1.140 at 
T ~ 2K. Above this temperature S(0) steadily decreases, and the renormaliza- 
tion is reduced to 1.098 at 4.2 K. The rather weak variation of the renormal- 
ization in this temperature interval means that the result for the renormalized 
temperature-dependent bulk susceptibility in the paramagnetic phase is close 
to the previous MF-result shown in Fig. 7 in Ref.l. In the presence of a mag- 
netic field, or in the ordered phase, the mean field is determined by (2.33). 
The change of this field, 5H, induced by the fluctuations, and the moments 
themselves are calculated fully self-consistently. Above Tc in the low-field limit 
5H/H MF = S(0). At higher fields, 5H goes through a maximum and then starts 
to decrease. At low temperatures the maximum value is about 0.3 kOe, and the 
maximum occurs when the moment is about half its saturation value (H MF ~ 4 
kOe). In the high field limit 5H/H MF vanishes, whereas 5H itself becomes small 
but not zero. It is uncertain whether this is due to the approximations made or 
not, but the non-zero value of 5H at infinite fields has no consequences. More 
importantly, the area determined by 5H as a function of gf^^{J x ) is found to 
agree quite accurately with the calculated energy change of the non-magnetic 
state induced by the fluctuations. In the low-temperature limit SU/N = —9.06 
fieV and (2.34) is satisfied to within 2-3%. Only in the high-temperature limit 
is the integral in (2.34) calculated to vanish somewhat faster than 5F(m = 0), 
indicating that the elastic contributions should be included in a more careful 
manner in this limit. 

At low fields the calculated renormalized magnetization is close to that 
predicted by the MF-model, corresponding to the behavior of the (zero-field) 
susceptibility. When the field is applied along the a-direction, the magneti- 
zation increases very rapidly and is close to its saturation value already at a 



12 



field of about lOkOe. At this field the l/z-renormalization is nearly quenched, 
corresponding effectively to an enhancement of the coupling constant. The cal- 
culated magnetization at low temperatures and intermediate fields is therefore 
increased in comparison with the result derived from the MF-model. At 4.2 K, 
when the field is applied along the c-axis, the (change of the) renormalization 
effect is small. In the other case, when the field is along the ct-axis at 1.6 K, 
the magnetization is calculated to be somewhat larger than predicted by the 
MF-model, and thus closer to the experimental behavior, as shown in Fig. 5. 
The figure only contains the results below 25kOe. The magnetization curves 
have been measured up to a field of 80 kOe. At the higher fields the results of 
the two models coincide and are in good agreement with experiment (see Fig. 6 
in Ref.l). The rapid quenching of the fluctuations produced by a field along the 
ct-axis at low temperatures is also illustrated by Fig.6, showing the field depen- 
dence of the excitations observed at (0 1) at 1.6 K. The width of the excitation 
band, which is close to the distance between the upper and lower excitation 
in this figure, is reduced by a factor of 5 at lOkOe, compared to its zero-field 
value. Above about 25kOe the excitation band is almost completely flat and 
the single-site fluctuations are of no importance. The change of the bare inter- 
action parameters means that the present model reproduces the experimental 
behavior somewhat more accurately than the previous model. 

Just above Tc, the entropy due to the electronic singlet-singlet states is 
vanishingly small, and the heat capacity in the ferrimagnetic phase is determined 
by the variation in the population of the nuclear levels as in the MF-case (thus 
the calculation of the heat capacity in the ordered phase does not rely on the 
relation (2.34)). The effects of the single-site fluctuations are included via the 
reduction of the mean field derived from (2.33). The result is compared with 
the experimental data in Fig.4 showing that the present model, in contrast 
to the MF-model, can almost account for the large jump in the heat capacity 
observed experimentally at Tc- The improved description of the heat capacity 
in the ordered phase is related directly to the fact that the ordered moment 
predicted in the zero-temperature limit is closer to the observed value. As 
shown in Fig. 7, the discrepancy is reduced by almost a factor of two, and the 
calculated moment at T = is now about 9% instead of 16% smaller than that 
determined 2 by the neutron-diffraction experiment. The larger value of the 
ordered moments implies an increase in the excitation energies and, as shown 
in Fig.8, the calculated excitation energies at 0.09 K are slightly above those 
derived from the MF-model. It is argued in Ref.l that the sample temperature 
at which the experimental results shown in Fig.8 were obtained might have been 
higher than indicated, because of the very large nuclear heat capacity at these 
temperatures. Here we find that we get a reasonable agreement between theory 
and experiment, if the sample temperature is assumed to be about 0.25 K. 
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IV. Discussion and conclusion 

The renormalization effects due to single-site fiuctuations have been included 
to leading order in the case of the singlet-singlet system H0F3. The expansion 
of the Matsubara Green function, applicable to this system, is considered to 
first order in l/z. In this order all the single sites are assumed to be placed in 
equivalent surroundings. The fluctuations in this effective medium, which de- 
rive self-consistently from the response of the single sites, affect the single-site 
Green function in a manner which may be determined by a cumulant expan- 
sion. An examination of the cumulants in the series of the single-site Green 
function shows that the usual 'unconditional' 1 /^-expansion could be improved 
substantially in a straightforward way. One additional advantage of the present 
procedure is that it is fully self-consistent, and the results therefore also behave 
properly close to a second-order phase transition. 

H0F3 is almost the ideal system for studying the effects of fluctuations. The 
electronic properties are nearly determined by the two lowest singlets alone at 
low temperatures, and the classical dipole-dipole coupling is a factor of ten 
larger than other dipole couplings in the system. The single-site properties 
are simple and the most important part of the two-ion interactions may be 
calculated directly. Furthermore, the long-range critical fluctuations close to the 
phase transition are expected to be of marginal importance, as they only lead 
to logarithmic corrections to the mean-field behavior. The effects of the l/z- 
fluctuations are calculated to be rather substantial in this system. For instance, 
the single-site susceptibility is found to be reduced effectively by up to 12% at 
the lowest temperature. In the comparison with the results of the MF-model 
in Ref.l, these differences are not always manifest, because the renormalization 
effects to some extent are included in this model via an effective adjustment 
of the interaction parameters. The renormalization effects vary slowly with 
temperature but may be quenched rather rapidly by applying a magnetic field, 
in which case differences may appear between the two models, since the bare 
interaction parameters are different. 

The comparison between theory and experiment has been improved system- 
atically by the inclusion of the effects due to the single-site fluctuations. The 
most striking improvement is found in the case of the heat capacity. Even so, 
the experimental heat capacity does not approach zero as closely as the the- 
oretical prediction at temperatures just above Tc, indicating that the critical 
fluctuations are of some importance, or alternatively that effects of the order 
(l/z) 2 should be included at the lowest temperatures. Although the theoretical 
magnetization curves, as a function of field or of temperature below Tc, are in 
better agreement with experiments than found in the MF-model, there are still 
some discrepancies. We do not expect that higher-order renormalization effects 
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can eliminate these deviations. An estimate of the bare electrostatic interaction 
between the quadrupole moments of the 4f-electrons on the different ions indi- 
cates that this coupling may be of some importance. However, this interaction 
is expected to be strongly shielded by the other electrons on the Ho- and the 
Cl-ions, and the good account of the heat capacity does not leave much room 
for any additional couplings. The discrepancies are systematic, but they are 
approaching a level where experimental uncertainties may be significant for the 
comparisons. Utilizing the correlation between the jump in the heat capacity 
at Tc and the zero-temperature moment, the comparison of experiment with 
theory suggests that this moment should be only a few per cent larger than 
calculated, between 5.2—5.4 /i B , to be compared with the neutron diffraction 
result 2 of 5.7 fx B ± 0.2 ji B . It is therefore apparent that a further refinement of 
the theoretical understanding of this unique magnetic system must depend on 
the performance of even more precise experiments. 
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FIGURE CAPTIONS 



Figure 1 The 4 Ho- and 12 F-ions in one unit cell of the orthorhombic structure of 
H0F3, projected on the a-c and the a-b plane. The arrows indicate the 
local x-axes lying in the a-c plane and making the angles ±9 with the 
a-axis, which are also the directions of the moments in the ordered phase. 
The figure shows the case where 9 = +25°. 

Figure 2 The final calculated density of states per Ho-ion in the low-temperature 
limit of the paramagnetic phase, J\f(e). The square- root singularities at 
the minimum and maximum energies are modified by the directional de- 
pendence of the long- wavelength Ai-mode (see Fig. 3). 

Figure 3 The dispersion relation of the singlet-singlet excitations along a* and c* 
in the paramagnetic phase of H0F3 at 1.6 K. The closed circles are the 
experimental results 1 obtained with the neutron-scattering vector along 
(h 0) and (0 0/), whereas the open circles are the results obtained along 
the energetically-equivalent directions {h 1) and (10/) (with h and / lying 
between and 1). The results are shown in a double zone representation, 
and the thin dashed lines indicate the Brillouin-zone boundaries. The 
solid lines are the theoretical predictions including the single-site fluctu- 
ations, and the thick dashed lines are the RPA-results derived from the 
MF- model in Ref.l. The labeling of the different modes, A\ - A4, close 
to the Bragg points is explained in Ref.l. The short solid lines marked 
A\ and A\ show the energies of the two ferromagnetic modes in the long 
wavelength limit, when the direction of the wave- vector is perpendicular 
to respectively (a 0) and (0 0/). In combination with the non-zero ex- 
perimental resolution width, the A\ mode in particular contributes much 
more strongly to the scattering cross-section than the A 1 mode close to 
(10 0). Resolution effects combined with the singular behavior of the long 
wavelength A r mode may also be important for explaining the difference 
between the theoretical and experimental energy of this mode. 

Figure 4 The low-temperature heat capacity of H0F3. The solid circles are the 
experimental results obtained by Bleaney et a/., 3 the dashed line is the 
magnetic contribution predicted by the MF-theory, 1 and the solid line is 
the result derived in the present approximation including the effect of the 
single-site fluctuations. Above 7-8 K the phonons start to contribute, as 
do the higher-lying 4f-levels on the Ho-ions. 
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Figure 5 The magnetic moment of H0F3 as a function of field along the a-axis at 
1.6 K and along the c-axis at 4.2 K. The experimental results of Bleaney 
et al. 3 have been corrected for the demagnetization field estimated from 
the shape of the samples. The mean-field predictions 1 are shown by the 
dashed lines, and the present calculations by the solid lines. 

Figure 6 The position of the inelastic scattering peaks at (0 01) as a function of 
magnetic field applied along the a-direction at 1.6 K. The circles are the 
experimental results 1 and the dashed and solid lines show the calculated 
energies of the A2, A3 and A3 modes, respectively, in the RPA 1 and when 
including the 1/z-renormalization. The A2 and A3 modes are near the 
minimum and maximum excitation energies, and the decreasing distance 
between the upper and lower lines indicates that the total excitation band- 
width rapidly declines as the field is increased. 

Figure 7 The solid lines show the calculated values of the angular moment, (J x ), 
and the nuclear spin, (I x ), of a Ho-ion as a function of temperature be- 
low Tc- The dashed lines are the corresponding results derived in the 
MF-case. 1 The filled circles are the experimental values of (J x ) deter- 
mined from the variation of the magnetic scattering intensities at (10 0) 
obtained by by Brown et at.? Their results have been scaled to agree 
with the magnetic moment of 5.7/ig, which they obtained from structure 
refinements at 70 mK. 

Figure 8 The dispersion relation along a* and c* in the ordered phase of H0F3 at 
90 mK. The meaning of the symbols is the same as in Fig. 3. The cross- 
hatched lines show the theoretical results at 0.25 K. 



18 



